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ABSTRACT 



Context. The Italian gamma-ray satellite AGILE has recently reported the detection of some variable high-energy 

sources likely of galactic origin. These sources do not have any obvious counterpart at lower energies. 

Aims. We propose that these sources are produced in proton-dominated jets of galactic microquasars. 

Methods. We develop a model for microquasar jets that takes into account both primary leptons and protons and all 

relevant radiative processes, including secondary particle emission and gamma-ray attenuation due to pair creation in 

the jet. 

Results. We obtain spectral energy distributions that correspond to what is observed by AGILE, with most of the power 
concentrated between 100 MeV and 10 GeV and reaching luminosities of 10 34-35 erg s _1 . We make detailed spectral 
predictions that can be tested by the Fermi gamma-ray telescope in the immediate future. 

Conclusions. We conclude that hadronic jets in galactic accreting sources can be responsible for the variable unidentified 
gamma-ray sources detected by AGILE. 
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1. Introduction 

The Italian gamma-ray satellite AGILE (Astro-rivelatore 
Gamma a Immagini Leggero) has detected several non- 
identified variable sources likely of galactic origin. These 
sources include the strong source AGL J202 1+4029 located 
in the Cygnus region, with the center of gravity of the error 
box at I = 78.01 deg and b = 2.19 deg (Longo et al. 2008), 
the variable source in the Musca region (error box centered 
at I = 312.2 deg and b = -0.3 deg, Pittori et al. 2008), and 
the high galactic latitude transient AGL J0229+2054, with 
the error box centered at I = 151.7 deg and b — —36.4 deg 
(Bullgarelli et al. 2008). The Cygnus source showed some 
significant re-brightening after its discovery (Giuliani ct al 
2008, Chen et al. 2008). Simultaneous X-ray observations 
with Super- AGILE did not show any counterpart in the 20- 
60 keV band. A steady and weak source was detected within 
the large error box by Swift/BAT (15-55 keV, see Ajello et 
al. 2008), but there is no clear relationship. Additional X- 
ray observations have been performed with XMM-Newton, 
without adding new clues (Pandel et al. 2008). Radio ob- 
servations with the Very Large Array (VLA) have shown 
no clear counterpart (Cheung 2008). 

Concerning the source in the Musca region, a poten- 
tial archival X-ray counterpart has been claimed on the 
basis of BeppoSAX observations dating from January 2001 
(Orlandini et al. 2008). The high-latitude transient, on the 
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other hand, might be a halo galactic source or a blazar (the 
AGN 1ES 0229+200 is at 61.3 arcmin from Gamma-Ray 
Imaging Detector error box centroid). 

The fact that these sources are highly variable im- 
plies that the high-energy radiation should be produced 
in a compact region. The absence of detection with Super- 
AGILE means that the ratio of gamma-ray to X-ray lu- 
minosities Lj/Lx should be >> 1. These characteristics 
recall those of the population of variable EGRET sources 
found in the Galactic Plane and in the Halo (Romero 
2001; Grenier 2001, 2004; Nolan et al. 2003). Actually the 
AGILE detections in the Cygnus region and the Musca 
region partially overlap with the location error box of 
the sources 3EG J2020+4017 and 3EG J1410-6147, respec- 
tively. It has been proposed that the unidentified variable 
gamma-ray sources sources at MeV-GeV energies might 
be high-mass microquasars with the emission dominated 
by inverse Compton up-scattering of UV stellar photons 
from the hot donor star (Kaufman-Bernado, Romero & 
Mirabel 2002, Bosch-Ramon, Romero & Paredes 2005). On 
the Galactic Plane, the donor star could be strongly ob- 
scured, rendering difficult its detection. These models, how- 
ever, predict a significant production of X-rays, something 
that is at odds with the new AGILE and Super-Agile ob- 
servations. Grenier et al. (2005) proposed that the variable 
high-latitude unidentified sources might be old, low-mass 
microquasars expelled long ago from the Galactic Plane 
or from globular clusters (see Mirabel et al. 2001). They 
also showed that external Compton models cannot account 
for the energetics required by the sources. Romero & Vila 
(2008) and Vila & Romero (2008) showed that "proton" 
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microquasars with low-mass donor stars might explain the 
halo EGRET sources through proton synchrotron radiation 
and photo-meson production. 

In this Letter we propose that microquasars with proton 
dominated jets can produce spectral energy distributions 
that satisfy all constraints imposed by AGILE observations 
and we make some predictions that can be used to test 
our proposal with the GLAST-Fermi satellite, and ground- 
based Cherenkov telescope arrays. 

2. The model 

The model we are going to present is based on the galactic 
jet model described by Romero & Vila (2008), but it in- 
corporates several refinements. The jet is assumed to carry 
most of the accretion power, in accordance with the dissi- 
pationless disk model of Bogovalov & Kelner (2005, 2008). 

A small fraction of the jet power is transformed into 
relativistic particles in a 'one-zone' acceleration region lo- 
cated close to the compact object (Khangulyan, Aharonian 
& Bosch- Ramon 2008). In this region, the magnetic field 
(obtained by equipartition) is extremely high (of the order 
of 10 7 G). This produces the immediate cooling of primary 
electrons and significant cooling of protons. In Fig. [T] we 
show the cooling rates for both primary electrons and pro- 
tons in the acceleration region, as well as the cooling and 
decaying rate of secondary muons and pions. The maximum 
energies of the primary particles are obtained equating the 
cooling rates to the acceleration rate, which is assumed to 
be produced by mildly relativistic shocks at the base of 
the jet through a first-order Fermi mechanism. In this work 
we adopt an acceleration efficiency of the order of 10 %, 
which rougly corresponds to a shock velocity of 0.3c and a 
mean free path similar to the Larmor radius. We see that 
electrons, even for such a high acceleration efficiency, reach 
only GeV energies, whereas protons can attain much higher 
energies, well into the PeV band. 

Our model takes into account: 1) synchrotron emission 
from both types (electrons and protons) of primary parti- 
cles, as well as emission from secondary leptons (electrons, 
positrons, and muons) and hadrons (charged pions), 2) in- 
verse Compton (IC) emission from all leptons in the strong 
total radiation field in the emitting region, 3) photo-pair 
and photo-meson production by both protons and pions, 
4) inelastic collisions between relativistic protons in the jet 
and the cold material that forms most of the same outflow 
(see Bosch- Ramon, Romero & Paredes 2006, for details of 
the physics of a cold matter dominated jet), 5) relativistic 
Bremsstrahlung from electrons and muons, and 6) internal 
photon absorption in the local radiation fields (calculated 
as in Aharonian, Khangulyan, & Costamante 2008), 7) re- 
injection of secondary pairs, which is negligible due to the 
low opacity of the proton-dominated jet (see Romero & Vila 
2008 for a discussion). The effect of losses of muons and pi- 
ons in a strong magnetic field on the resulting gamma-ray 
and secondary lepton spectra has been recently discussed 
by Reynoso & Romero (2008) and we adopt here their treat- 
ment. For the other processes we follow the classical formu- 
lae (Blumenthal & Gould 1970), the photo-meson produc- 
tion treatment already used by Romero & Vila (2008), and 
the expressions given by Kelner, Aharonian, & Bugayov 
(2006) for pp interactions (this formalism is only valid for 
proton energies above 0.1 TeV; at lower energies we use the 
standard treatment). 



The ratio of relativistic protons to electrons a in the 
jet is unknown and is used as the basic free parameter in 
our model. We have calculated a number of models for dif- 
ferent acceleration efficiencies, and a set of standard pa- 
rameters: black hole mass ~ 10 Mq, accretion rate ~ 10 37 
erg s _1 (i.e., ~ 10~ 2 of the Eddington luminosity), injec- 
tion cx 25 — 1-5 , viewing angle 30°, jet bulk Lorentz factor 
r = 1.5, and location of the acceleration zone at zq rj 10 s 
cm (see Romero & Vila 2008). 

The transport equation is solved for all types of parti- 
cles. At steady state it reads: 
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where t csc is the particle escape time from the acceleration 
region (t csc « Az/vj et ) and Q is the injection function that 
can be normalized in accordance with the energy budget of 
relativistic particles through: 



L e . B = I d 3 r 



dE e _p E e ^p Q e,p(-^e,p) ■ (2) 



The minimum kinetic energy is taken to be of the or- 
der of the rest mass energy of the corresponding particle. 
The equation for pions and muons has an additional term 
N(E, z)/ta cc , that takes into account the decay of the par- 
ticles on a timescale idee that depends on the energy in the 
lab system. 

In general, variability can be obtained introducing a 
variable injection Q(E,t). Once the particle distributions 
are known, the radiative output can be calculated as men- 
tioned above. The internal attenuation results from photon 
annihilation: 7 + 7 — > e~ + e + . The opacity for a gamma 
ray of energy E 7 is: 



r(E 7 ) = - 



1 Jtth 



(1 — it) cr(E 7 , e, it) 



xn(e, z) du de dl, (3) 



where n(e, z) is the photon number density at energy e and 
location z, it = cost?, 1? is the angle between the momenta 
of the colliding photons, I is the 7-ray photon path, and the 
cross section for the interaction is given by (e.g. Levinson 
2006): 



a^(E y , e, 0) = — a T (l-p 2 ) 



(3-/3 4 )ln(i±|)-2/3(2~ / 3 2 ) 



(4) 



Here, (3 is the speed of the electron/positron in the center 
of momentum frame, i.e., 



(l-/3 2 ) = 



2m2 c 4 



< p < 1. 



(1 - U )£ 7 e' 

The threshold energy eth is defined by (3 = 1 with ■& = 
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Fig. 1. Acceleration and cooling rates at the base of the jet for primary protons and electrons, and secondary pions 
and muons, calculated for representative values of the model parameters (proton-to-lepton energy ratio a = 1000, and 
primary injection spectral index a = 1.5). The value of the magnetic field is B ~ 10 G, and the acceleration efficiency 
parameter T) is indicated. 
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Fig. 2. Spectral energy distributions of a proton-dominated microquasar (a = 1000). Each panel corresponds to a different 
acceleration efficiency (77 = 0.1 on the left, r\ = 0.01 on the right). Fermi point source sensitivity is indicated by a thick 
blue line. 



3. Results total power of 10 37 erg s _1 . We show the different contribu- 

tions from all significant cooling processes for all primary 
In Fig. [2] we show the spectral energy distribution (SED) and secondary particles, in the case of two different ac- 
calculatcd for a proton dominated (a = 1000) jet, with a 
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celeration efficiencies 77 (0.1 and 0.01, from left to right, 
respectively). Bremsstrahlung and IC losses are negligible, 
since synchrotron radiation is the dominat cooling channel 
for leptons. The different acceleration efficiencies not only 
affect the maximum energy of the radiation but modify the 
whole shape of the SED, since the dominant synchrotron 
losses follow the square of the energy and depend signif- 
icantly on the mass of the particles, resulting in a global 
photon redistribution. In both cases, however, the peak of 
the SED is determined by proton synchrotron radiation, 
followed by pion-synchrotron emission. For the higher ef- 
ficiency the synchrotron peak is sharper, reaching almost 
10 35 erg s _1 . In the case of a lower efficiency, the peak 
is slightly above 10 34 erg s _1 . In the first case most of the 
emission is concentrated in the range 10 8 — 10 12 eV, whereas 
in the second it is between 10 8 and 10 10 eV, with a soft slope 
beyond 10 9 eV. In both cases there is a high ratio L 7 /Lx, in 
accordance with what is inferred from AGILE observations. 
Soft X-rays, due mainly to electron synchrotron radiation, 
are at the level of 10 32 erg s _1 . The hard X-ray component 
is dominated by muon synchrotron emission. Contrary to 
models with equipartition (a = 1) in relativistic particles, 
photo-meson production is not significant in proton domi- 
nated jets, since the synchrotron field is relatively weak (see 
Romero & Vila 2008). 

The quoted luminosities correspond to the flux reported 
by AGILE and the upper X-ray limits, at a distance of the 
order of ~ 2 kpc. And accurate determination of the dis- 
tance can be used to constrain the energy budget of the 
sources. At TeV energies the luminosities (due mainly to 
pp interactions) are at the level of 10 32 — 10 33 erg s _1 , and 
hence the sources are not detectable by current Cherenkov 
atmospheric telescopes, although they could fall within 
the sensitivity of upgraded instruments like HESS II and 
MAGIC II. 

4. Discussion and summary 

A crucial feature of a synchrotron-proton dominated jet is 
that a very strong magnetic field is necessary to produce 
detectable radiation in 7-rays. This means that the accel- 
eration region should be located very close to the compact 
object (at ~ 10 8 cm in our models). Intrinsic absorption is 
then not very important since the low energy fields, which 
are responsible for the opacity to gamma-ray propagation, 
are weak. In models with a high content of primary elec- 
trons these effects, at the base of the jet, are very signifi- 
cant leading to a complete suppression of all emission above 
~ 100 GeV (Romero & Vila 2008). Such models produce a 
huge amount of X-rays, something that is not observed in 
the unidentified MeV-GeV sources. 

In our model we introduce a hard injection spectrum 
in order to achieve a strong contrast between leptonic 
and hadronic peaks. A softer injection would reduce the 
L 1 /Lx ratio. We notice that the losses in the high mag- 
netic field strongly affect the overall leptonic particle spec- 
trum. Observations with the LAT instrument of the Fermi 
satellite will allow us to determine the photon spectrum of 
these sources in the range 100 MeV — 100 GeV, allowing a 
better determination of the physical parameters. 

We remark that IC models are possible for these type of 
sources only if the radiation is produced in a region with a 
significantly smaller (~ 10 — 100 G) magnetic field, see for 
example Punsly ct al. (2000). Such a field would occur at 



distances of 10 12 — 10 13 cm from the compact object. How 
the injector of relativistic particles could be located there 
will be discussed elsewhere (see, nonetheless, Bosch-Ramon 
2007). 

In summary, we propose a complete lepto/hadronic jet 
model to explain the unidentified variable AGILE sources. 
This model assumes a strong component of relativistic pri- 
mary protons and takes into account all radiative processes 
that might occur at the base of the jets. The predicted 
SEDs are in accordance with what we know about of these 
sources. The jet model is independent of the nature of the 
donor star, so it could explain both low- and high-latitude 
galactic sources. Fermi observations will allow us to de- 
termine better constraints on the spectral features, then 
making possible inferences about the actual conditions in 
the sources. Development of models for 7-ray production 
far from the compact object and models that take into ac- 
count external effects is in progress. 
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